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Abstract 

We compute vacuum expectation values of products of fermion 
^ ■ bilinears for two-dimensional Quantum Chromodynamics at finite 

flavored fermion densities. We introduce the chemical potential 
as an external charge distribution within the path-integral ap- 
proach and carefully analyse the contribution of different topo- 
logical sectors to fermion correlators. We show the existence of 
chiral condensates exhibiting an oscillatory inhomogeneous be- 
havior as a function of a chemical potential matrix. This result 
is exact and goes in the same direction as the behavior found in 
QCD^ within the large ./V approximation. 
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1 Introduction 



In order to understand the structure of the QCD vacuum one should anal- 
yse possible mechanisms for chiral symmetry breaking and the formation of 
fermion condensates. The existence of such correlators can be understood 
as the result of condensation of pairs of particles and holes and it can have 
interesting implicancies in particle physics and cosmology. For example, a 
color nonsinglet condensate may be related to superfluidity and color su- 
perconductivity of cold quark matter at high fermion densities 0. In this 
respect the results of Deryagin, Grigoriev and Rubakov || are of particular 
importance. Analysing the large N c limit of QCD these authors have shown 
that the order parameter for chiral symmetry, the quark condensate (i/ji/j), 
is at high quark densities inhomogeneous and anisotropic so that, regarding 
the order parameter, the ground state of quark matter has the structure of 
a standing wave. 

Two-dimensional models like the Schwinger model and QCD 2 provide a 
natural laboratory to test these phenomena since, although simplified, the 
basic aspects (chiral symmetry features, non-trivial topological sectors, etc) 
are still present and exact calculations can be in many cases performed. 

An analysis of two-dimensional QED at finite density was originally pre- 
sented in 0-0. More recently, studies on this theory showed that 
inhomogeneous chiral condensates do exist as a result of the contribution of 
non-trivial topological sectors. 

Extending our work on QED2 @ we analyse in the present paper vacuum 
expectation values of products of local bilinears ijj(x)ip(x), at finite density 
for two-dimensional Quantum Chromodynamics with flavor. Using a path- 
integral approach which is very appropriate to handle non-Abelian gauge 
theories, we show that the multipoint chiral condensates exhibit an oscilla- 
tory inhomogenous behavior depending on a chemical potential matrix. Our 
results are exact and, remarkably, go in the same direction as those revealed 
in four dimensions using the 1/N C approximation to QCD 0]. 

To study the effect of finite fermion density in QCD 2 a chemical poten- 
tial may be introduced. Within the path-integral approach this ammounts to 
consider a classical background charge distribution in addition to that pro- 
duced by topologically non-trivial gauge configurations. Concerning this last 
point, it is well-known that in two space-time dimensions the role of instan- 
tons is played by vortices. In the Abelian case, these vortices are identified 
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with the Nielsen-Olesen solutions of the spontaneously broken Abelian Higgs 
model ||. Also in the non- Abelian case, regular solutions with topological 
charge exist when symmetry breaking is appropriately achieved via Higgs 
fields 0-[|lO]]. In both cases the associated fermion zero modes have been 
found |J-!J. 

Properties of the vortex solutions and the corresponding Dirac equation 
zero-modes are summarized in section 2. We then describe in sections 3 
and 4 how topological effects can be taken into account within the path- 
integral formulation leading to a compact form for the partition function in 
the presence of a chemical potential. Our approach, following ref. WM, starts 



by decomposing a given gauge field belonging to the n th topological sector in 
the form 

A^x) = Ajjfi + A? + a, (1) 

Here Afi> is a (classical) fixed gauge field configuration belonging to the n th 
class, Af^ is the background charge field taking account of the chemical 
potential, and a M is the path-integral variable which represents quantum 
fluctuations. Both A e * 1 and a M belong to the trivial topological sector and 
can be then decoupled by a chiral rotation with the sole evaluation of a 



Fujikawa jacobian ||15|| . This last calculation can be easily performed since it 
is to be done in the trivial topological sector. 

The complete calculation leading to the minimal non-trivial correlation 
functions of fermion bilinears is first presented for multiflavour QED 2 (Sec- 
tion 3) and then extended to multiflavour QCD 2 (Section 4). In both cases 
the oscillatory behavior of correlators as a function of the chemical poten- 
tial is computed, the result showing a striking resemblance with the QCD^ 
answer obtained within the large N c approximation |§. We summarize our 
results and conclusions in section 5. 



2 Zero Modes 

Topological gauge field configurations and the corresponding zero-modes of 
the Dirac equation play a central role in calculations involving fermion com- 
posites. We sumarize in this section the main properties of vortices, the rele- 
vant topological objects in the model we shall consider, both for the Abelian 
and non-Abelian cases. We also present the corresponding Dirac operator 
zero-modes. 
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2.1 The Abelian case 



In two-dimensional Euclidean space-time, topologically non-trivial configu- 
rations are available since Nielsen and Olesen || presented their static z- 
independent vortex. In the U(l) case the topological charge for such a con- 
figuration, working in an arbitrary compact surface (like a sphere or a torus) 
is defined as 

l.Jd 2 xe, u F^ = neZ (2) 

A representative gauge field configuration carrying topological charge n 
can be written as 

A^=ne^A(\x\) (3) 

with v4(|x|) a function which can be calculated numerically (an exact solution 
exists under certain conditions on coupling constants, |§). The adequate 
boundary conditions are 

A(0) = , lim A(\x\) = -1 (4) 

\x\— >oo 

There are \n\ zero-modes associated with the Dirac operator in the back- 
ground of an A™ configuration in a suitable compactified space-time Jl4| 



(For the non-compact case see JTTJ] ) . For n > (n < 0) they correspond 
to right-handed (left-handed) solutions (t]l) which in terms of light-cone 
variables z = xo + ix\ and z = xo — ix\ can be written in the form 

Vn = ( Zmh{ Q ZrZ) ) (5) 

V ™ = ( z' m h-\z,z) ) (6) 
where m — 0, 1, . . . , \n\ — 1, 

h(z,z)=exp[<f ) ( ~ n \\z\)] (7) 

and 

d ^»)(|*|) = nA(|*|). (8) 



d\z 
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2.2 The non-Abelian case 



As in the Abelian case, two-dimensional gauge field configurations car- 
rying a topological charge n e Z N can be found for the SU(N) case. As 
explained in ref.[]TB[ the relevant homotopy group is in this case and not 
Z as in the U(l) case. 

Calling ip the angle characterizing the direction at infinity, a mapping 
g n {}p) G SU(N) belonging to the n th homotopy class (n = 0, 1, . . . , N — 1) 
satisfies, when one turns around a close contour, 

g n (2n) = exp(^)g n (0) (9) 

Such a behavior can be achieved just by taking g n in the Cartan subgroup 
of the gauge group. For example, in the SU(2) case one can take 

g n (<p) =exp[Va i (y?)] (10) 

with 

^„(2tt) - fi n (0) = 2n(2k + n) (11) 

Here n — 0, 1 labels the topological charge and fc £ Z is a second integer 
which conects the topological charge with the vortex magnetic flux (Only for 
abelian vortices both quantities coincide). 

We can then write a gauge field configuration belonging to the n th topo- 
logical sector in the form 

4"> = iA{\x\) g-%g n (12) 
with the boundary conditions 

A(0) = , lim A(\x\) = -1 (13) 

\x\— >oo 

These and more general vortex configurations have been thoroughfully stud- 
ied in 0-[T§. 

Concerning zero-modes of the Dirac operator in the background of non- 
Abelian vortices, they have been analysed in refs.[T2|]-[[nj. The outcome 
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is that for topological charge n > (n < 0) there are Nn (N\n\) square- 
integrable zero modes t)l (t]r) analogous to those arising in the Abelian case. 
Indeed, one has 



^ = {z' m hl\zrz)) (15) 

with 

h{z,z) =exp[^ n \\z\)M] (16) 

and 

M = ldiag(l,l,...,l-iV) (17) 

Here i,j = 1,2, ... ,N and m = 0, 1, . . . , \n\ — 1. The pair (m,i) labels 
the jV|n| different zero modes while j corresponds to a color index. Due to 
the ansatz discussed in refs.||lU||-[]r6] for the non- Abelian vortex, the function 



^(|z|) appearing in eq.fllBD coincides with that arising in eqs. ([?])-(§) for the 
abelian vortex. 

As it happens in the abelian case, the partition function of two dimen- 
sional Quantum Chromodynamics only picks a contribution from the trivial 
sector because det^fA^]) = for n ^ (see eq. ([77]) below). In contrast, 
various correlation functions become non-trivial precisely for u ^ thanks 
to the "absortion" of zero-mode contributions when Grassman integration is 
performed. 

It is our aim to see how these non-trivial correlators are modified when 
a fermion finite density constraint is introduced, comparing the results with 
those of the unconstrained (zero chemical potential) case. As explained in 
the introduction, we are motivated by the results of Deryagin, Grigoriev and 
Rubakov || in four dimensional QCD. They were able to show, in the large 
N c and high fermion density limits, the existence of oscillatory condensates 
(the frequency given by the chemical potential) which are spatially inhomo- 
geneous. For QED 2 the same oscillatory behavior was found approximately 
in |P and confirmed analytically in [^], by examining an arbitrary number 
of fermion bilinears for which the exact /i-dependence of fermionic correla- 
tors was computed. In order to improve our understanding of the large N c 
results found in QCD4, we shall extend in what follows our two-dimensional 
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approach to the non-Abelian case but before, we shall consider the case of 
flavored QED 2 as a clarifying step towards multiflavor QCD 2 . 



3 Multiflavour QED 2 

We developed in ref . a path-integral method to compute fermion compos- 
ites for Abelian gauge theories including chemical potential effects. In this 
section we briefly describe our approach while extending our treatment so as 
to include flavour. We then leave for section 4 the analysis of the non-Abelian 
multiflavour QCD 2 model at finite density. 



(i) Handling the chemical potential in the Abelian case 

We start from the Lagrangian 

L = -^F flv F flv + $(i@+4-iMjo)'4> (18) 

where ip is the fermion field isospinor. A chemical potential term has been 
included by considering the diagonal matrix M. defined as 

M = diag(//i . . . n Nf ) (19) 

where Nf is the total number of flavors and /!& are Lagrange multipliers 
carrying a flavour index, so that each &;-fermion number is independently 
conserved. The corresponding partition function is defined as 

Z\pti . . . fi Nf ] = J V^VipVA^expi- J d 2 x L). (20) 

Since our interest is the computation of fermionic correlators, we have to 
carefully treat non-trivial topological configurations of the gauge fields which 
have been seen to be crucial in the obtention of non-vanishing condensates, 
see refs. [P^] -P9[. Then, following the approach of refs.|i^-[|T9|j, we decompose 



gauge field configurations belonging to the n th topological sector in the form 

A,(x) = A^(x) + a,(x) (21) 

where is a fixed classical configuration carrying all the topological charge 
n, and a M , the path integral variable, accounts for the quantum "fluctuations" 
and belongs to the trivial sector n — 0. 
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As it is well-known |[20|| , the chemical potential term can be represented 
by a vector field A e * 1 describing an external charge density acting on the 
quantum system. Indeed, taking A e * 1 as i times the chemical potential matrix 
(see eqs. (|19D and (|22|)) it corresponds to a uniform charge background for 



each fermionic flavor. As explained in 0, it is convenient to first consider a 
finite length (21) box and then take the I — > oo limit. In this way translation 
symmetry breaking associated to the chemical potential becomes apparent 
and simultaneously, ambiguities in the definition of the finite density theory 
are avoided. When necessary, we shall follow this prescription (see ref . Q for 
a discussion on this issue). We start by defining 



so that the Dirac operator 

can be compactly written as 
with 



A e f = -iM 5„ , (22) 

i@+ 4 - iM-fo (23) 

ip+ 4' (24) 

1 ext 



A'^ = A, + Af (25) 

We shall now proceed to a decoupling of fermions from the chemical 
potential and the a M fluctuations following the steps described in M for the 
case of only one flavor. In that case, we wrote 



a„ = -e 



and made a chiral rotation to decouple both the <p — f] fields together with 
the chemical potential. In order to include Nf flavors in the analysis, one 
has to replace (0, 77) — > (4>,r])lf and n — ► M. as we shall see below. Then, 
we can straightforwardly apply what we have learnt for one flavor [0 in the 
multiflavor case. The change of variables accounting for the decoupling of 
fermions from the field together with the chemical potential is given by 

■0 = exp[7 5 (cj)(x)lf + iMxi) + vq(x)lf] x 

(27) 

ip — x exp[7 5 ((j)(x)lf + iMxi) - irj(x)lf] 
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^ = (-Z0C7) t/- 1 (28) 

where 

U = exp[ 75 (01/ + iMxx) + ^l f ] (29) 

For notation compacteness we have included in /& the external field A e * 1 
describing the chemical potential term. From here on we choose the Lorentz 
gauge to work in (which in our notation corresponds to r\ — 0). 

After transformation (^7|) the resulting Dirac operator takes the form 

%jp = ift+ 4 {n) + i -> i0+4 (n) . (30) 

The jacobian associated with the chiral rotation of the fermion variables can 
be easily seen to be |7| 

J = exp J d 2 x (01/ + iMx^U^ + 20< n >)) (31) 
where (p^ is defined by 



Together with eq. (|27|) we consider the change in the gauge- field variables a M 
so that 

T>a„ = A FP 5(ri)V<i)Vri (32) 

with App = det □ 

As thoroughly analysed by Actor |20| , A e * 1 does not correspond to a pure 



gauge. Were it not so, the introduction of a chemical potential would not 
have physical consequences and this would be the case in any space-time 
dimensions. In fact, one cannot gauge away L chem by means of a bounded 
gauge transformation. As explained in j7|, the chiral rotation which decouples 
the chemical potential, although unbounded can be properly handled by 
putting the system in a spatial box, then introducing adequate counterterms 
and finally taking the infinite volume limit. 

After the decoupling transformation, the partition function, can be writ- 
ten in the form 

Z = NY. VxD X V<P exp(-£<#) (33) 

n J 

where is the effective action in the n th topological sector, 
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S% = I d 2 x m + 4 {n) )x-^ 2 J d'x^f + e^F^) 

-0 / d 2 x(F$f j d 2 x (01 / + ^x 1 )D(0 + 20(")) + S c (34) 

The usual divergency associated to the electromagnetic energy carried by 
fermions has to be eliminated by a counterterm S c ||. I n °ur approach the 
divergency manifests through the term iA4x\0(j)^ in eq. . Putting the 
model in a finite length box and appropriately adjusting S c yields to a finite 
answer. The counterterm is the Lagrangian counterpart of the one usually 
employed in the Hamiltonian approach to handle this problem |4j]. In the 
canonical formulation of QFT this is equivalent to a redefinition of creation 
and annihilation operators which amounts to a shift in the scale used to 
measure excitations. 

As we have mentioned, a fermionic chemical potential ammounts to intro- 
ducing a finite external charge (i.e. at the spatial boundaries) into the theory. 
In these conditions, it can be proved that massless QED 2 (and QCD 2 ) at 
finite density remains in the Higgs phase. To show this one may compute the 
string tension following for instance the procedure described in ref. [£TJ|: one 
starts by integrating out the fermion fields (or, equivalently, the bosons in the 
bosonized version) in order to derive the effective action for the gauge fields. 
One can then compute the Wilson loop to calculate the energy of a couple 
of (static) external charges for a theory containing also dynamical 'quarks'. 
Now, since zero modes kill the contributions of non-trivial topological sectors 
to the partition function, screening can be tested using the effective action in 
the trivial topological sector. In fact, one can see that for vanishing fermion 
masses the string tension vanishes. In order to discuss these issues in multi- 
flavour QED 2 at finite density, let us note that after integration of fermions 
in eq. fl34|) , the resulting effective action for the gauge field can be written as 

H 

/ I | %J\A. r 

d2x ^2 F l + ~ tr f^ J d " X ^0, ( 35 ) 

where S c has cancelled the divergent term, as explained above. Then, at this 
stage there is no divergency to deal with and we can perform our calculation 
in the whole Euclidean space. Choosing the Coulomb gauge a\ = 0, ap- 
propriate to derive the static potential between external charges, we obtain 
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from (0), (after integration by parts using zero boundary conditions) the 
following effective Lagrangian 

£eff = ^(^o) 2 + 7^0 - tr /7^ a ( 36 ) 

In order to analize the force between charges, let us pass to Minkowski space- 
time, making ao — > iao so that the corresponding effective Lagrangian in 
Minkowski space reads 



£m = -^(^i«o) 2 - -^ajj + tr f ^a (31 ) 



To determine the electrostatic potential between two external charges ±e', 
we may couple to the gauge field the proper charge density 

p^ 1 ) = e '{5{x l + I) - Six 1 - I)) (38) 

so that the complete effective Lagrangian becomes 

£ = C M - pa (39) 

The resulting equation of motion takes the form 

e 2 e 2 

dla a + —tr } M = p (40) 

7T Z7T 

and its solution reads 
e' 

a (x ) = - — (exp(— mix 1 + l\) — exp(— mjx 1 — l\)) + trjM. (41) 



where m = e^Nj/n. The energy of the two test charges at a distance 21 
each other, is given by 

V(l) = - f dx 1 p(x 1 )a (a; 1 ) (42) 



2 

and we obtain, at finite fermion density, the usual screening potential 

V(l) = —(l-exp(-2ml)) (43) 
10 



with no modification due to the presence of the chemical potential whose 
contribution trivially cancels. We then conclude that in massless multiflavour 
QED 2 at finite density, any fractional charge e' is screened by integer massless 
charges. 

To get a deeper insight into these results, let us note that, as it is well- 
known, only the trivial topological sector contributes to the partition function 
for massless fermions (the contribution of non-trivial sectors being killed by 
zero-modes). Then, deriving the potential between external charges as we 
did above or computing the Wilson loop W as well, one finds that screening is 
not affected by the presence of the chemical potential term (the Wilson loop 
calculation yields W = 1 pip. One could argue that, as it happened with the 
two test charges, the external charge background associated to the chemical 
potential term is itself also screened by massless dynamical fermions. After 
all it is only in the properties of fermion condensates that topological sectors 
enter into play and it is through its contributions that the chemical potential 
manifests itself. One can understand this issue as follows: The topological 
structure of the theory is determined at the boundaries (recall that in order 
to calculate the topological charge one just uses § A^dx^ = 2nn) and the 
corresponding n-charge configurations are responsible for the non-triviallity 
of the correlators but not affecting the partition function. It is then precisely 
when computing condensates, that the charges at the boundaries associated 
with the chemical potential manifest. 

To rephrase this analysis, note that with the choice of the counterterm dis- 
cussed above, the effective action written in terms of the decoupled fermions 
does not depend on the chemical potentials /z^. This does not mean that 
this term has no physical consequences. In fact, M. reappears when com- 
puting correlation functions of fermion fields, once ip and ip are written in 
terms of the decoupled fields \ an d X through eq.(|55|). We shall see in the 
following sections how fermionic correlators are changed, exhibiting oscilla- 
tory inhomogeneities in the spatial axes which depend on M.. The fact that 
zero modes make certain v.e.v.'s not to vanish, leads to a highly non-trivial 
dependence on the chemical potentials. 

(ii) The Correlation Functions 

The introduction of a flavor index implies additional degrees of freedom which 
result in Nf independent fermionic field variables. Consequently, the grow- 
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ing number of Grassman (numeric) differentials calls for additional Fourier 
coeficients in the integrand. 

It is well known that each coeficient is related to the quantum numbers 
of the chosen basis, which is normally builded up from the eigenf unctions 
of the Dirac operator. As we have for one flavor, one has that n of these 
eigenfunctions are zero-modes, implying a vanishing fermionic exponential. 
Hence, in order to make Grassman integrals non-trivial, one has to insert 
several bilinears depending on the number of zero modes. When the path- 
integral measure contains Nf independent fermionic fields instead of one, 
the number of composite insertions is multiplied by Nf in order to saturate 
Grassman integration algebra, with some selection rules which will become 
apparent below. 

For the sake of brevity let us readily give the result for general correlation 
functions of p points with arbitrary right and left insertions 



N f r k s k 




(44) 



k=l i=l j=V 



where 



4K) = viMviM , 



(45) 



k=l 




d 2 x 0D(0 + (n) )] x 



k=l i=l 
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I] j^W-W exp[- / d 2 x x k ^+ 4 {n) )x k ] (46) 

where w\ is the space component of u>\ We see from eq. fl46|) that the chemical 
potential contribution is, as expected, completely factorized. Concerning the 
bosonic integral, it can be written as 

/Nf r k Sfe 
d 2 x <p^U<p^} exp[-2^(^0 (n) (w'') - E <p{w j ))] 

k=l i=l j=V 

N f Pk Pk' 

xexp[-2 £ EX^i^VV)] (47) 

k,k'=l i=l i=l 

with 

_1 (w 4 , w- 7 ) = fT ('H' u;4 - ^ J |) + hi(c|u/ - w j \). 
The fermionic path-integral determines the topological sectors contributing 
to equation fl4~4l) . More precisely, once the correlator to be computed has 
been chosen, Grassman integration leads to a non-zero answer only when 
the number of right insertions minus the number of left insertions is the 
same in every flavor sector. It means that — = t V7c, where t is 
the only topological flux number surviving the leading sumatory in eq. (|46|) . 
(Notice that mixed flavor indices in the elementary bilinear are avoided, 
i.e. we are not including flavor- violating vertices, in accordance with QED4 
interactions). It is important to stress that each term explicitly including the 
classical configuration of the flux sector cancels out. Consequently, classical 
configurations only appear through by means of their global (topological) 
properties, namely, through the difference in the number of right and left 
handed bilinears fT4}| . 

To conclude, we give the final result for the general correlator defined in 
eq. fl4"4j ) making use of the explicit form of Abelian zero modes 

N f r k s k 7 

<nn4(^n^)> = (-^) p 

k=l i=l j=l> 47F 

N f r k s k N f Pk Pk' 

exp[2iJ2^k(J2 w l~J2 w i)] II exp[-4^X>* e i ln ( c K- w J '|)] 

k=l 8=1 j=l' k>k'=l i=l j=l 

Nf p k p k i 

exp[- E E E e l e 3 K (m\w t - w j \)} (48) 

k,k' i=l j=l 
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(see Ref.0 and [|24| for details). 

In order to clearly see the meaning of this expression, let us show the 
result for the simplest non-trivial flavored correlation functions including 
mixed right and left handed insertions 

n 

2cos[//i(zi - x 1 - y x ) - ^ 2 wi}{s\(x)s l + (y)s\(z)s 2 + (w)) 1 + 
2cos[//i(yi - X! - zt) - fi2W 1 ](s 1 + (x)s 1 _(y)s 1 + (z)s 2 + (w)) 1 + 
2cos[/ii(xi -z\ -yi) -// 2 wi](sUz)s+(y)s+(z)s+M):i., (49) 

n 

2cos[/ii(xi - yi) - // 2 (zi - u;i)](s+(x)s^(y)s 2 _(2;)s^(w)) + 
2cos[/ii(xi - yi) + // 2 (^i - ^i)](s 1 f (x)si(y)s+(2;)s^(u7))o + 
2cos[/ii(xi +yi) +At 2 (^i + w 1 ))(s 1 + (x)s 1 _(y)s 2 _(z)s 2 + (w)} 2 (50) 

These expresions make apparent: (i) How the topological structure of the 
theory exhibits itself through the existence of non-trivial vacuum expectation 
values of fermionic bilinears. (Notice that those on the right hand side are the 
only surviving terms of the whole sumatory). (ii) In the multiflavor case, the 
path-integrals are non-zero only when the number of right insertions minus 
the number of left insertions are identical in every flavor sector, (iii) The 
sum over spatial coordinates dramatically exhibits the translation symmetry 
breaking discussed above, (iv) The fixing of various fermion densities implies 
a somehow reacher spatial inhomogeneity of the results with respect to the 
one flavor case that we have analyzed in M , in the sense that now the "angles" 
depend on various chemical potentials, (v) Another difference with respect 
to the one flavor case, concerns the trivial cancellation of logarithms coming 
from bosonic and fermionic integration respectively. Now, this cancellation 
occuring for one flavor, does not take place anymore, see eq.([48|). 

4 Multiflavour QCD 2 

In the present section we consider two dimensional SU (N c ) Yang-Mills gauge 
fields coupled to massless Dirac fermions in the fundamental representation. 
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Due to the non-Abelian character of the gauge symmetry, gluons are charged 
fields that preserve color flux at each vertex. Since a colored quark density 
is not a quantity to be kept constant, no chemical potential related to color 
should be considered but only that associated with the global symmetry that 
yields fermion number conservation. Hence, we first include one chemical 
potential term and then consider a different lagrange multiplier for each 
fermionic flavor. 

Let us stress that once the topological effects arising from vortices are 
taken into account and the chemical potential behavior of fermion correlators 
is identified, we do not pursue calculations in the bosonic sector (neither we 
consider the inclusion of Higgs scalars, necessary at the classical level for the 
existence of regular vortex solutions). As we shall see, the boson contribution 
to the fermion condensate just factorizes and all the chemical potential effects 
can be controlled by calculations just performed within the fermionic sector. 

(i) Handling the Chemical Potential in QCD2 
We start from the massless QCD 2 (Euclidean) Lagrangian 

l = nid^' + A^ a tf'i, - + ^Kv F %- ( 51 ) 

where we have included a chemical potential term in the form 

Lchem = -i/JLIp^lp (52) 

in order to take care of the fermion density constraint. Here a = 1 . . . N% — 1, 
and q — 1 . . . N c . The partition function reads 

Z\fi]= J Vrj/D^VA^expl- J d 2 x expL]. (53) 

Again, one can decouple the chemical potential by performing an appropriate 
chiral rotation for the fermion variables. Indeed, under the transformation 

V> = exp(?/i7 5 :ri)x 

i> = X exp(i/i7 5 a;i) (54) 
the fermion Lagrangian becomes 
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so that the chemical potential completely disappears from the fermion La- 
grangian. As we have seen, chiral transformations may generate a Fujikawa 
jacobian which has to be computed using some regularization procedure. For 
example, using the heat-kernel regularization one introduces a resolution of 
the identity of the form 

1 = lim exp(- P(a) 2 /M 2 ). (56) 

M— >oo 

where D^a) (a G (0, 1)) is an interpolating Dirac operator such that -D M (a = 
0) = D„[A, fi] and D^{a = 1) = D^[A}. 

After some standard calculation |G3| one ends with a Jacobian of the form 



J = exp (^ie^/An J d 2 x da tr c [fiXiF flu (a)] S j (57) 

where tr c is the trace with respect to color indices and 

F^(a) = F« u (a)t a , a = 1, 2, . . . , N 2 C - 1 (58) 

Now, the color trace in eq.(^) vanishes and then the chiral Jacobian is in 
fact trivial, 

J = 1 (59) 

We can then write the partition function (^) after the fermion rotation 
defined in eq.(|27|) in the form 

Z\pi\ = J VA^VxVxeM- J d 2 xL) (60) 

As we have seen in the Abelian case, although fx is absent from the r.h.s. of 
eq . (|60|) one should not conclude that physics is independent of the chemical 
potential. For correlation functions of composite operators which are not chi- 
ral invariant, the chemical potential will reappear when rotating the fermion 
variables in the fermionic bilinears. As in the Abelian case, this happens 
when computing v.e.v.'s of products i]j(x)iJj(x) 



(ii) Correlation functions in QCD2 with chemical poten- 
tial 

Our main interest is the computation of fermionic correlators containing 
products of local bilinears ipi/j(x) for which non-trivial topological gauge field 
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configurations, and the associated Dirac operator zero-modes, will be crucial 
to the obtention of non- vanishing results as explained in refs.0,[14]-[19|. 



As in section 3, we start by writing a gauge field belonging to the n th 
topological sector, in the form 

A#x) = Af*\x) + dfa) (61) 

where A a ^ is a fixed classical configuration (as described in section 2.2) 
carrying all the topological charge n, and a", will be the actual integration 
variable which belongs to the trivial sector n = 0. Then, we decouple the 
a M field from the fermions through an appropriate rotation (the calculation 
of the Fujikawa Jacobian being standard since the decoupling corresponds 
to the topologically trivial sector). Now, it will be convenient to choose the 
background so that 

A a + {n) = (62) 
In this way, the Dirac operator takes the formQ 



and we are left with the determinant of this operator once fermions are 
integrated out 

Z[/A = £/ Va^M^F'AA^ + a]] det p[A™ + a}. (64) 

As before, we have introduced a sum over different topological sectors. Now, 
we shall factor out the determinant in the classical background so as to 
control the zero mode problem. Let us start by introducing group valued 
fields to represent A^ and a M 

a + = iu~ d + u (65) 

a_ = idiyvd-V^ 1 )^ 1 (66) 
A<l ] = idd^d- 1 . (67) 



1 We are using 70 = o\ and 71 = — 02. 
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Consider first the light-cone like gauge choice |23 

A_ = A W (68) 

implying 

v = I. (69) 
In this gauge the Dirac operator flBUD reads 

W (n) + a]| Ic = ( a _ + A (n) o ) ( 7 °) 



where subscript Zc means that we have used the gauge condition (|68|). One 
can easily see (for example by rotating the + sector with w" 1 while leaving 
the — sector unchanged) that 

det p[A {n) + a] |fc = A/" det x exp(W[u, A {n) ]). (71) 

Here W[u, A- n '\ is the gauged Wess-Zumino-Witten action which in this case 
takes the form 

W[u, A (n) ] = + — tr c f (Px^d+u^dd-dr 1 ) (72) 

47T J 

and W[iz] is the Wess-Zumino-Witten action 

W[u] = — tr c / d 2 xdu,u~ 1 du,u H tr c / d 3 y (■u~ 1 9j-u)('U~ 1 9j-u)(-u~ 1 (9fcM). 

2n J ' An Jb 

(73) 

Note that in writing the fermion determinant in the form (fn]), the zero- 
mode problem has been circumscribed to the classical background fermion 
determinant. 

One can easily extend the result ( fTTD to an arbitrary gauge, in terms of 
the group- valued fields u and v defined by eqs.(|65D-(|6^), by repeated use of 
the Polyakov-Wiegmann identity [p6 



W[pq] = W[p] + W[q] + — tr c [ d 2 x{p~ l d + p) {qd^q' 1 ) (74) 

An J 

The answer is 

det I/)[A {n) +a]= A/" det l/)[A {n) } x exp(S eff [u, v; A {n) }) (75) 
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S eff [u,v;A^} = W[u,A {n) } + W{v] + ^-trJd 2 x{u- 1 d + u)d(vd^v- 1 )d- 1 

4:71 J 

+—tr c [d 2 x(d- 1 d + d)(vd„v- 1 ). (76) 

47X J 

Once one has the determinant in the form ([75Q, one can work with any gauge 
fixing condition. The gauge choice (^) is in principle not safe since the 
corresponding Faddeev- Popov determinant is A = det D a 3 \A^\ implying 
the possibility of new zero-modes. A more appropriate choice would be for 
example A + = 0, having a trivial FP determinant. In any case one ends with 
a partition function showing the following structure 

Z = £det(#[A< n >]) fva„A5(F[a}) 

exp a,} - ± J d?xFl v [A {n \ «j) (77) 

Concerning the divergency associated to the external charge distribution, we 
have learnt from the Abelian case that one has to carefully handle this term 
in order to define excitations with respect to the external background. In 
section 3 we have seen that it came from the interaction of A ext with Fj$, 
appearing in the fermionic jacobian. Performing a similar calculation in the 
present case we would find the non-Abelian analogue of this term with tr c 
acting on it. As we have mentioned above, this color trace operation implies 
the vanishing of the corresponding divergency so that no counterterm might 
be added in QCD 2 , meaning that the relevant vacuum is properly defined. 

As we have seen, the Lagrangian for QCD 2 at finite density can be writ- 
ten in terms of //-rotated fields which hide the chemical potential from the 
partition function. This result however, does not exhausts the physics of the 
theory in the sense that correlation functions do depend on \i. Actually, it 
will be shown that the chemical potential dependence appears as a factor 
multiplying the result for correlators of the unconstrained theory. For this 
reason, we shall first describe the computation of vacuum expectation values 
of fermion bilinears in the /t = case and then consider how this result is 
modified at finite fermion density. Hence, we proceed with the analysis of 
v.e.v's of products of bilinears like XX- Let us start by noting that with the 
choice (B^) for the classical field configuration, the Dirac equation takes the 
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form 



where ( is defined as 

X+ = dvd~ x ( + 

X- = u-\- (79) 
so that the Lagrangian in the n th flux sector can be written as 

L =xP[a + AW]x = CJd+C- + C P-[A {n) K + 

= CD[A^}(. (80) 

In terms of these new fields, the bilinears XX take the form 

Xx = Cudvd- 1 ^ + C+dv^d^u- 1 ^- (81) 

We observe that the jacobian associated to flT9|) is nothing else but the ef- 
fective action defined in the previous section by eq. (|76|) . Hence, an explicit 
expression for the non-Abelian correlators reads 

(xxix 1 ) ■ --XXix 1 )) = E / Va » A 6 ( F M exp[-S of f(A^\a)) 

n J 

B q ^(x x ) . . . B qm (x l ) C-X+ix 1 ) . . . CC+(x l ) + B^{x x ) . . . 
B -^iPi(x l ) CTC+ (x 1 ) ■ ■ .C+'C-(x l ) + B^{x l ) . . . 
B- lqi -^(x l - l )B- lqiPl (x l ) C-X+ix 1 ) . . . C+"- 1 C- , " 1 (a:'~ 1 )C+ , C- 
+ ... (82) 

where the group-valued field B is given by 

B = udvd' 1 . 

For brevity we have written the gauge field measure in terms of the orig- 
inal fields although for actual calculations in the bosonic sector one has 
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to work using u and v variables and proceed to a definite gauge fixing. That 
is, the measure should be written according to 



Da^ — ► VuT>vJb{u, v, d) 

and then the gauge condition and Faddeev-Popov determinant should be 
included (For example, in the light-cone gauge a + = 0, u = 1 and the FP 
determinant is trivial). Finally, notice that we have obtained a general and 
completely decoupled result, from which one sees that due to color degrees 
of freedom, the simple product that one finds in the Abelian case becomes 
here an involved summatory. 

Now that we have an expression for correlators in the unconstrained case, 
let us include the chemical potential in our results. Recall that in this theory 
the partition function is (see eq.([60|)) 



Z = J VA^VxV X exp (- J d 2 x X (i@+ 4)x + ^F^ u 



(83) 



where x, x represent the fermion fields after the chiral rotation (54) which 
eliminated the chemical potential from the Lagrangian. Since fermionic bi- 
linears can be written as 

one has 

(ipip) = exp(2i/j,x 1 ) (x+x+) + exp(-2?>ri) (x-X-)- (84) 

It can be easily seen that the same factorization occurs when flavor is 
introduced. The corresponding transformation for the fermion field isospinor 
is now 

$ = x exp(0/U c 7 5 :ri) (85) 

and the bilinear v.e.v takes in this case the form 

(ipip) = exp(2iMl c x 1 )(x+X+) + exp(-2zA^l c x 1 )(x-X-)- (86) 

We shall then include from here on flavor degrees of freedom with the corre- 
sponding constraint on each fermion density. Since in this case one deals with 
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Nf fermions coupled to the gauge field, we can use the fermionic jacobian 
we have computed for one flavor to the power Nf while the bosonic mea- 
sure remains untouched. In the light-cone gauge it can be easily seen that 
the effective bosonic sector now involves N c — 1 massive scalars, their mass 
depending on flavor and color numbers by means of a factor (2N C + Nf) 1 / 2 
with respect to the abelian counterpart (There is also the same number of 
unphysical massless particles p7p. 



As we have previously explained, the Dirac operator has \n\N c zero modes 
in the n th topological sector, this implying that more fermion bilinears are 
needed in order to obtain a non-zero fermionic path-integral. Moreover, since 
the flavor index implies a factor Nf on the number of Grassman coeficients, 
the minimal non-zero product of fermion bilinears in the n th sector requires 
\n\N c Nf insertions. 

Since the properties of the topological configurations are dictated by those 
of the torus of SU (N c ), one can easily extend the results already obtained for 
QED2. In particular, the chirality of the zero modes is dictated by the same 
index theorem found in the Abelian theory, this implying that in sector n > 
(n < 0) every zero mode has positive (negative) chirality. In this way, the 
right (left) chiral projections of the minimal non-zero fermionic correlators 
can be easily computed. One gets 

Nf N <= |n| 1 r m 

(nnn^W))" = w vuVvj * e-w^ 

k q i " JGF 



N f No M Nc , . r-~ , s \ 

nnn e ^r^f ) ( / vm j^ w k m 



where 



B' 9 k ' pik {x) = exp(2i/i fc x 1 ) u Piq {x){dvdr l ) qk {x), (88) 

C, + = and D[A^ n ^] stands for the Dirac operator in the r.h.s of eq. (|80|) . 
We have used the notation for the partition function since it is com- 



pletely determined within the n = sector, see eq.(|77j). We have showed 
every color and flavor indices explicitly indicating sum and product opera- 
tions. The GF label stands for the gauge fixing. The action S^hJu, v, d) = 
NfSwzw( u , v i d) + SMaxweii(u, v, d) is given by the full gluon field A^ n >(d) + 
a(u,v), and yields a high order Skyrme-type lagrangian W7 . 
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Let us consider N c = 2 and Nf = 2 in order to present the simplest 
illustration for the last expression. The minimal fermionic correlator then 
looks 



n 

i N c =2 2 

7^ E I[exp[2ifi k (x 1 + Vif] VuVvJ E 

L p,q,r,sk=l JGF 

Bl' pq (x k )B 2 / s (y k ) J V( k V( k er<^ Kk ftktftf) (y *). (89) 

The fermionic path-integral can be easily done, resulting in the product 
of eigenfunctions discussed in the sections above, as follows 

/ VCuDCk eS^ A(1) ^ (f(t(x k ) C/CAy k ) = det/pL4«]) x 

+# 2 >'Sf 2)9,fc (^ fc )# 1)r 'U° ,1)S ' fc (2/ fe ))- (90) 

Here det /(.DL^ 1 )]) is the determinat of the Dirac operator defined in eq.flSTp 
omitting zero-modes and (e.g.) r)( 0,1 ^ q ' k (x k ) is a non-Abelian zero-mode as 
defined in section 2, with an additional flavor index k. Concerning the bosonic 
sector, the presence of the F^ u (Maxwell) term crucially changes the effective 
dynamics with respect to that of a pure Wess-Zumino model. One then 
has to perform approximate calculations to compute the bosonic factor, for 
example, linearizing the U transformation, see ||17|| . In any case, once this 
task is achieved for the \i = model, the modified (finite density) result can 
be obtained in an exact way. 



5 Summary 

We have presented the correlation functions of fermion bilinears in multi- 
flavour QED 2 and QCD 2 at finite fermion density, using a path-integral 
approach which is particularly appropriate to identify the contributions aris- 
ing from different topological sectors. Analysing correlation functions for 
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an arbitrary number of fermionic bilinears, we have been able to determine 
exactly its dependence with the chemical potentials associated to different 
flavor indices. As stressed in the introduction, our work was prompted by 
recent results by Deryagin, Grigoriev and Rubakov || showing that in the 
large N c limit, condensates of QCD in four dimensions are inhomogeneous 
and anisotropic at high fermion density. 

Two-dimensional models are a favorite laboratory to test phenomena 
which are expected to happen in QCD4. In fact, an oscillatory inhomo- 
geneous behavior in (ipip) was found in the Schwinger model [fj using oper- 
ator bosonization and then the analysis was completed by finding the exact 
behavior of fermion bilinear correlators in [^]. Here we have extended this 
analysis in order to include flavor and color degrees of freedom within a path- 
integral scheme which makes apparent how topological effects give rise to the 
non-triviality of correlators. 

Remarkably, the oscillatory behavior related to the chemical potential 
that we have found with no approximation, coincides exactly with that de- 
scribed in || for QCD 4 within the large N c approximation (appart from 
the anisotropy that of course cannot be tested in one spatial dimension). 
In particular, the structure of the multipoint correlation functions, given by 
eqs. ([48D and (|S7D, shows a non-trivial dependence on spatial coordinates. 
This makes apparent that the ground state has, at finite density, an involved 
structure which is a superposition of standing waves with respect to the order 
parameter. Being our model two-dimensional, we were able to control the 
chemical potential matrix behavior in an exact way so that we can discard 
the possibility that the formation of the standing wave is a byproduct of 
some approximation. This should be considered when analysing the results 
of ref.|| in d — 4 dimensions, where one could argue that use of a ladder 
approximation as well as the fact of neglecting effects subleading in 1/N C 
play an important role in obtaining such a behavior. 

Several interesting issues are open for further investigation using our ap- 
proach. One can in particular study in a very simple way the behavior of 
condensates at finite temperature. The chiral anomaly is independent of 
temperature and plays a central role in the behavior of condensates through 
its connection with the index theorem. Therefore, one should expect that 
formulae like ( ^8|) or (|87| ) are valid also for T > 0. Of course, v.e.v.'s at 
li = in the r.h.s. of this equation, should be replaced by those computed at 
finite temperature and hence the issue of zero-modes in a toroidal manifold 
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should be carefully examined (see e.g. |24|]). At the light of recent results 



concerning QCD2 with adjoint fermions JH, ^8fl-||29|| it should be of interest 
to extend our calculation so as to consider adjoint multiplets of fermions. 

Finally, it should be worthwhile to consider massive fermions and compute 
fermion correlation functions at finite density, via a perturbative expansion 
in the fermion mass following the approach of P0| . We hope to report on 
these problems in a future work. 
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